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NOTATION 


A  Amplitude  of  angular  variation  oscillation 

D  Propeller  diameter 

Ct  Propeller  blade  chord  length  at  r/R=0.7 

J  Advance  ratio  =  U^^/nD 

1,  Downstream  propeller  wake  surface  contact  length 

L  Lift  force  on  blade  /  span  of  blade 

n  Propeller  rps 

r  Radial  coordinate 

R  Propeller  radius 

Re  Reynolds  number  =  UD/v 

Re7  Reynolds  number  based  on  0.7R  blade  chord  =  C7  [U2^  +(0.7nD7i)2]OJ  /  v 

Carriage  velocity 

x  Streamwise  distance  from  plane  of  propeller  blade  midchords 

Zfc*  Depth  of  hub  axis  from  free  surface 

P  Wake  spreading  angle 

r  Vortex  strength  ■  L/pU 

v  Kinematic  viscosity 

0,  Angular  displacement  of  propeller 

0,  Angular  displacement  of  motor  drive  shaft 

p  Density 

Angular  velocity  of  propeller 
<o.  Angular  velocity  of  motor  drive  shaft 
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ABSTRACT 


Propeller  wake  flow  visualization  was  carried  out  for  several  depths  of 
submergence  below  a  free  surface.  The  wake  was  visualized  by  laser  light 
sheet  illumination  of  the  blade  tip  vortex  helix  seeded  with  fluorescein  dye. 
Effects  of  the  free  surface  on  the  propeller  wake  were  seal  at  the  shallower 
submergence  depths.  Instability  in  the  helix  was  apparent  at  much  shorter 
distances  downstream  than  has  previously  been  observed  in  cavitation 
experiments.  The  dominant  cause  of  this  instability  is  determined  to  be 
caused  by  the  use  of  inclined  universal  joints  in  the  propeller  drive  system. 

A  simple  test  modification  is  described  for  the  elimination  of  this  influence. 
Estimation  of  Reynolds  number,  submergence  depth,  and  propeller  loading 
influences  on  the  propeller  wake  are  given. 


ADMINISTRATIVE  INFORMATION 

The  work  described  in  this  report  is  part  of  die  Surface  Ship  Wake  Consortium 
sponsored  by  the  Office  of  Naval  Research  (ONR),  under  Program  Element  61153N,  Task 
Area  BT02301N1,  and  performed  under  die  David  Taylor  Research  Center  (DTRC)  work 
unit  1-1504-200  (FY  88).  Additional  funds  were  provided  by  the  6.2  Ship  and  Submarine 
Technology  Program,  Program  element  62543N,  Task  area  SF  43421  for  analysis  of  the 
experimental  results  and  publication  of  this  report. 

INTRODUCTION 

Recent  studies  suggest  that  the  interaction  of  turbulence  in  a  ship’s  wake  with  the 
free  surface  may  define  characteristics  of  some  parts  of  the  ship  hydrodynamic  signature 
(Lyden,  et  al.,  1985).  As  a  major  producer  of  turbulence  in  the  ship’s  wake,  the  propeller 
deserves  careful  examination  to  determine  both  the  character  of  its  turbulent  wake  and  the 
resulting  influence  this  wake  has  on  the  free  surface.  Velocity  measurements  made  at 
David  Taylor  Research  Center  (DTRC)  have  shown  a  lack  of  Made  rate  periodicity  in  the 
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propeller  wake  beyond  five  propeller  diameters  downstream  (Blanton  and  Fish,  1988).  This 
degradation  of  blade  rate  dependance  suggests  either  turbulent  diffusion  of  the  individual 
blade  wakes  or  an  instability  in  flow  geometry  disturbing  the  wake  axisymmetry. 
Photographs  of  propellers  cavitating  in  water  tunnels  have  shown  very  symmetric  helical 
structure  in  the  propeller  wake.  Water  tunnel  photographs  are  typically  limited  in 
downstream  viewing  range  to  x/D=2.0  (Kuiper,  1979).  In  addition,  the  tunnel  walls  may 
provide  some  stabilizing  effect  on  the  propeller  wake  geometry.  In  order  to  examine  the 
interaction  of  the  propeller  wake  with  a  flee  surface,  the  structure  of  the  wake  should  be 
identified  and  understood.  This  understanding  is  of  fundamental  importance  in  validating 
numerical  prediction  codes  aimed  at  simulating  the  propeller  wake  interaction  with  the  flee 
surface.  Due  to  the  complicated  nature  of  this  flow  field,  a  visual  perspective  on  the  wake 
evolution  was  deemed  necessary  prior  to  any  detailed  wake  survey  measurements. 

This  report  reviews  flow  visualization  data  collected  on  a  propeller  wake  operating 
at  two  depths  below  a  flee  surface.  High  contrast  black  and  white  video  coverage  showing 
laser  light  sheet  illumination  of  the  tip  vortex  helical  wake  were  recorded  for  a  variety  of 
operating  conditions.  Screen  images  are  included  here  with  associated  comments.  A 
preliminary  report.  Fish  and  Blanton  (1988),  provided  initial  presentation  and  analysis  of 
this  experiment. 


EXPERIMENTAL  ARRANGEMENT 

The  David  Taylor  Research  Center  140  foot  Towing  Basin  was  used  for  the 
experiments  presented  in  this  report.  Although  the  carriage  war  operated  at  speeds  from  1 
to  4  ft/sec  (0.305  to  1.22  m/sec),  the  data  presented  here,  except  where  noted,  were  taken 
at  a  carriage  speed  of  1  ft/sec  (0.305  m/sec)  for  maximum  dye  clarity. 
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A  four  bladed,  12  in.  diameter  propeller  (DTRC  #3563)  was  used.  This  propeller 
was  modified  by  milling  a  1/8  in.  (0.3175  cm)  wide  slot  in  each  blade  face.  The  milled 
slot  was  then  covered  with  tape  to  form  a  passage  for  dye  from  the  hub  to  the  Made  tip  as 
shown  in  Fig.  1.  The  0.002  in  (0.1  mm)  thickness  of  die  tape  was  assumed  small  enough 
to  maintain  the  original  blade  face  geometry  and  resulting  propeller  performance. 

A  tri-strut  and  pod  configuration  was  designed  and  built  to  minimize  both  the  free 
surface  disturbance  and  the  propeller  inflow  asymmetry  while  maintaining  structural  rigidity 
(see  Fig.  2).  The  propeller  drive  motor  was  remotely  located  on  the  carriage  and 
connected  to  the  hollow  prop  shaft  through  a  3/8  in.  (0.95  cm)  angle  drive  shaft  positioned 
just  forward  of  the  center  strut.  U-joints  on  either  end  of  the  angle  drive  shaft  were 
aligned  to  minimize  angular  velocity  variations  (no  visual  differences  in  angular  velocity 
were  initially  apparent;  see  WAKE  HELIX  INSTABILITY  section). 

A  variable  speed  positive  displacement  pump  supplied  a  fluorescein  dye  solution  to 
the  pod.  The  fluorescein  dye  was  illuminated  using  a  laser  light  sheet,  positioned  as 
shown  in  Fig.  3.  The  light  sheet  was  created  using  a  110  mW  argon-ion  laser  directed 
into  a  30  Hz  vibrating  mirror.  This  sheet  was  fixed  (not  moving  with  the  carriage)  in  the 
vertical  plane  and  aligned  with  the  propeller  axis. 

A  Sony  3/4  in.,  30  frame  per  second  video  tape  recorder  was  used  to  record  images 
from  a  high  sensitivity  Mack  and  white  camera  mounted  on  the  basin  side  as  shown  in 
Fig.  3.  The  images  shown  here  are  35  mm  photographs  of  the  video  screen. 
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DYE  INJECTION  FLOW  VISUALIZATION 


The  dye  injection  rate  was  examined  to  determine  its  influence  on  die  formation  of 
the  tip  vortex.  Care  was  taken  to  balance  the  need  for  maximum  dye  flow  rate  producing 
extended  wake  visualization  time  with  the  need  for  minimum  dye  flow  rate  disturbance  of 
the  flow  in  the  blade  tip  region.  Selection  of  an  "optimum"  dye  delivery  flow  rate  was 
made  after  recording  wake  images  under  a  variety  of  dye  flow  rates  and  propeller  rotation 
speeds.  "Optimum"  here  was  defined  as  the  minimum  flow  rate  giving  sufficient  dye  in 
the  vortices  for  sharp  contrast  to  one  diameter  downstream.  Dye  pump  speed  variation 
about  this  small  flow  rate  had  no  apparent  influence  on  the  vortex  wake.  The  resulting 
flow  rate  of  0.394  in’/sec  (6.43  ml/sec)  gave  an  estimated  radial  velocity  on  the  blade  face 
exit  point  of  16.8  in  ./sec  (43  cm/sec).  This  exit  velocity  is  approximately  one  third  of  die 
circumferential  velocity  at  the  tip  and  is  oriented  in  the  same  direction  as  the  redial  flow 
induced  by  the  tip  vortex.  Recall  that  the  dye  slot  is  on  the  face  or  "pressure”  side  of  the 
blade. 

Dye  injection  rates  were  increased  only  in  select  high  speed/low  J  cases  in  order  to 
enhance  the  video  contrast.  In  these  cases,  no  apparent  changes  were  noticed  in  the  flow, 
only  sharper  contrast  in  the  images.  The  large  variation  in  exit  velocity  of  the  dye  caused 
by  the  flow  rate  change  indicated  that  under  these  circumstances,  the  momentum  transfer 
from  the  dye  exit  point  had  a  negligible  effect  on  the  formation  of  the  tip  vortex. 

Care  must  be  taken  in  interpreting  dye  roll-up  as  discrete  vortices.  As  shown  by 
Hama  (1962)  and  recently  by  Cimbala  et  al.  (1988),  streakline  dye  patterns  show  the 
integrated  motion  of  the  marked  fluid  from  its  point  of  introduction.  In  the  current 
experiment,  however,  local  rates  of  change  of  the  dye  roll-up  were  sufficient  to  indicate 
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persistent  vortical  structure.  The  continuous  change  in  the  dye  roll-up  provided  an  accurate 
marker  of  the  fluid  velocity.  Comparison  with  the  velocity  field  induced  by  discrete  vortex 
lines  leads  to  accurate  location  estimates  of  the  vorticity  centers. 

WAKE  HELIX  INSTABILITY 

The  most  striking  feature  illustrated  in  the  video  coverage  in  this  experiment  was 
the  twisting  of  vortex  structures  in  pairs  resembling  a  2 -dimensional  vortex  "pairing" 
phenomenon  (Ho  and  Huerre,  1984)  in  the  plane  of  the  laser  light  sheet.  This  flow 
phenomenon  will  be  referred  to  as  vortex  twisting  and  will  be  described  further  in  die 
following  paragraphs.  Fig.  4  shows  a  time  sequence  of  the  helix  in  0.1  second  increments. 
The  propeller  loading  corresponds  to  J=0.5.  The  dye  spots  marking  the  vortex  cores  can 
be  seen  to  pair  up  and  rotate  about  each  other  as  the  propeller  moves  further  to  the  right 
of  the  photographs.  If  the  wake  were  to  remain  stable  and  symmetric,  the  dye  spots  would 
remain  roughly  equidistant  from  each  other,  appearing  similar  to  a  discretised  vortex  sheet. 
The  movement  of  these  dye  spots  therefore  indicates  a  loss  of  symmetry  caused  by  an 
unsteady  influence  on  the  propeller  wake.  This  type  of  vortex  motion  appeared  in  all  of 
the  operating  conditions  of  this  experiment  except  at  the  non-loaded  advance  ratio.  Details 
of  any  differences  are  described  in  later  sections  of  this  report.  Although  previous  analytic 
work  indicated  die  inherent  instability  of  helical  vortices  (Widnall  1972),  propeller  tip 
vortex  cavitation  observations  in  water  tunnels  have  shown  stable  helix  forms.  To 
determine  die  nature  of  the  observed  vortex  twisting,  an  evaluation  of  possible  influences 
on  propeller  wake  instability  was  conducted. 
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Sources  of  unsteady  perturbation  on  propeller  wake  flow  fields  are  listed  below: 

a.  unsteady  propeller  inflow 

b.  non-uniform  spatial  propeller  inflow 

c.  flow  asymmetry  caused  by  die  free  surface 

d.  unsteady  rotation  of  propeller 

e.  blade-to-blade  variations  in  loading 

It  should  be  noted  that  item  d  can  be  modeled  as  a  combination  of  items  a  and  b. 

Item  e,  blade-to-blade  loading  variations,  could  be  caused  by  geometric  variations 
between  the  blades,  or  by  variations  in  the  boundary  layer  flow  pattern  over  similarly 
shaped  blades.  Dimensional  measurements  of  the  propeller  however,  showed  very  uniform 
blade  geometries.  In  addition,  no  roughness  variations  were  found  from  blade-to-blade 
which  might  cause  boundary  layer  transition/separation  variations. 

Since  these  experiments  were  conducted  in  a  still  water  towing  tank,  the  turbulence 
level  was  observed  to  be  extremely  low  and  spatially  uniform  prior  to  each  run.  The 
influence  of  die  strut  wakes  on  the  wake  helix  could  not  be  neglected,  but  appeared  to  be 
insignificant  judging  by  the  similarity  of  vortex  twisting  in  die  lower  region  near  the 
propeller  (no  strut  effect)  to  that  of  the  upper  region  near  the  propeller  (See  Fig.  4).  The 
effect  of  strut  wakes,  if  any,  would  be  seen  in  this  region.  Any  inflow  unsteadiness  was 
therefore  suspected  to  be  due  to  carriage  motions  or  propeller  drive  unsteadiness.  Although 
vibrations  on  the  carriage  were  not  measured,  the  same  form  of  unsteady  flow  was 
observed  over  wide  variations  in  both  carriage  speed  and  blade  rate  frequencies.  The 
similarity  of  the  flow  unsteadiness  under  these  conditions  implied  a  negligibly  low 
influence  of  unsteady  carriage  motions. 
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To  determine  if  an  unsteady  perturbation  of  the  flow  field  was  being  created  by  the 
propeller  drive  system,  careful  measurements  of  the  angular  motion  of  the  propeller  were 
made.  Fig.  5  shows  results  of  a  comparative  measurement  of  the  angular  positions  of  die 
motor  drive  shaft  0,  and  the  propeller  9,.  The  two-cycle  nature  of  the  curve  is  a 
characteristic  of  the  non-constant  angular  transmission  of  the  two  inclined  universal  joints. 
Despite  their  ease  in  design  application,  universal  joints  in  this  inclined  fashion  were  found 
to  cause  an  undesirable  unsteady  angular  velocity. 

An  analysis  shows  how  this  small  unsteady  angular  velocity  variation  can  dominate 
the  downstream  interaction  of  the  propeller  tip  vortices  and  could  give  rise  to  the  vortex 
twisting  observed  in  this  experiment.  The  key  to  the  analysis  is  the  coincidence  of  the 
four  lobes  (upper  and  lower  peaks)  in  the  angular  velocity,  to,  with  the  positions  of  the 
four  blades  of  the  propeller.  A  typical  sequence  of  blade  passages  through  a  fixed  angular 
position  can  help  in  visualizing  the  significance  of  this  result. 

a.  Assume  that  blade  #1  passes  through  the  vertical  position  when  the  propeller 
is  operating  at  the  positive  peak  of  03  variation. 

b.  Blade  #2  will  pass  through  the  vertical  position  90  degrees  later,  but  now  the 
propeller  is  operating  at  the  negative  peak  of  angular  velocity. 

c.  Blade  #3  passes  vertical  180  degrees  later  and  the  propeller  is  again  at  the 
positive  peak  of  CD  similar  to  when  blade  #1  passed  through  the  vertical 
position. 

d.  Finally  blade  #4  passes  through  vertical  with  the  minimum  to  similar  to  when 
blade  #  2  passed  through  vertical. 
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In  other  words,  the  180  degree  separation  of  every  other  blade  synchronizes  with  the  180 
degree  separation  in  positive  (or  negative)  peaks  in  the  angular  velocity  curve.  The  net 
result  is  an  oscillation  (at  each  fixed  angular  position)  of  the  instantaneous  value  of 
propeller  advance  ratio.  In  this  experiment,  the  light  sheet  focuses  the  viewer’s  attention 
on  two  discrete  positions:  top  and  bottom.  The  visible  vortices  shed  into  the  wake  will 
therefore  alternate  in  strength  and  position  based  on  the  advance  ratio,  J,  of  die  propeller  at 
the  instant  the  blade  passes  through  the  sheet.  With  these  conditions  defined,  one  can 
derive  a  vortex  interaction  scenario  for  comparison  with  the  observations. 

The  fluctuations  in  J  caused  by  the  <D  oscillation  give  rise  to  changes  in  the  axial 
compression  of  the  vortex  helix  geometry.  A  tip  vortex  will  experience  an  imbalance  in 
induced  velocity  from  its  neighboring  blade  tip  vortices  due  to  their  unequal  spacing.  If 
the  closest  neighbor  vortex  is  downstream,  the  net  induced  motion  will  be  radially  inward. 
The  significance  of  this  inward  motion  arises  from  the  high  radial  gradient  in  axial  velocity 
in  the  tip  vortex  region.  This  gradient  cruises  filaments  of  lower  radial  position  to  be 
convected  downstream  faster  than  their  partners  at  higher  radial  position.  This  convection 
can  be  visualized  in  both  the  3-Dimensional  images  of  Fig.  6  and  die  laser  light  sheet 
images  (see  Fig.  4).  Note  in  Fig.  4  the  relative  movement  of  the  dye  spots  labeled  "1" 
and  "2".  Despite  their  similar  motion  to  vortex  "pairing",  it  should  be  remembered  that  the 
twisting  process  as  shown  in  Fig.  6  is  very  3-dimensional. 

The  continued  interaction  of  the  vortices  once  they  have  been  shed  was  seen  to  be 
dependant  on  the  propeller  loading  and  is  discussed  further  in  the  PROPELLER  LOADING 
VARIATION  section. 
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The  dominance  of  unsteady  angular  velocity  in  the  propeller  wake  will  be 
eliminated  in  future  experiments  through  die  use  of  constant  velocity  joints.  Despite  the 
inclusion  of  this  unsteady  co  variation,  the  results  obtained  from  the  present  study  will  be 
useful  in  die  study  of  propeller  wakes  near  a  free  surface. 

PRESENTATION  OF  RESULTS 
REYNOLDS  NUMBER  DEPENDENCE 

Experimental  measurements  such  as  (Cook,  1972)  have  shown  that  tip  vortices  can 
be  modelled  with  finite  radius  viscous  cores,  outside  of  which  the  flow  behaves  in  an 
inviscid  manner.  The  existence  of  a  viscous  core  leads  one  to  examine  die  Reynolds 
number  dependence  of  results  utilizing  visualization  of  these  vortices.  A  series  of  cases 
was  therefore  recorded  for  several  carriage  speeds  while  varying  the  rpm  to  maintain  the 
same  advance  coefficient,  J.  The  resulting  increase  in  flow  speed  over  the  blades  could 
therefore  be  examined  without  changing  the  flow  angle  of  attack.  This  allowed  observation 
of  wake  variations  due  to  viscous  effects.  This  comparison  was  made  only  with  the  deep 
model  submergence  (minimizing  Froude  number  effects  due  to  flee  surface  waves). 

Fig.  7  isolates  the  upper  row  of  wake  vortices  at  J=0.8  for  three  different  Re^ 
ranging  from  «1.3  x  10s  to  *2.6  x  105.  Although  the  rate  of  dye  diffusion  is  greater  for 
die  higher  Re^,  the  trajectory  and  pairing  of  the  vortices  are  essentially  the  same.  This 
indicates  that  under  the  conditions  of  this  test,  the  Reynolds  number  effects  were 
overshadowed  by  the  effect  of  the  angular  velocity  variation  described  previously  in  die 
WAKE  HELIX  INSTABILITY  section. 

Reynolds  numbers  (Re^)  for  the  propeller  varied  from  *5  x  10*  to  *1.5  x  10s  for  a 
carriage  velocity  of  1  ft/s.  Meyne  (1972)  reports,  for  similar  4-bladed  propellers,  the 
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pressure  side  remained  laminar  up  to  Re^  =  8  x  10s,  and  the  suction  side  remained 
laminar  between  Re&7  =  2  x  105  and  Re^  =  3  x  10s,  depending  on  advance  coefficient. 
This  suggests  that  die  visualization  data  for  the  1  ft/sec  carriage  speed  corresponds  to 
laminar  flow  over  the  entire  propeller. 

As  the  Reynolds  number  is  increased,  the  laminar  tip  vortices  will  transition  to 
turbulent  flow  at  smaller  axial  distances  downstream  (see  Fig.  7).  In  order  to  reveal  any 
effects  of  laminar  to  turbulent  vortex  transition  on  die  wake  interaction  process,  however, 
the  propeller  drive  modifications  to  remove  unsteady  0)  variations  will  be  needed.  The 
constant  propeller  angular  velocity  should  give  rise  to  a  more  symmetric  vortex  wake, 
increasing  the  importance  of  viscous  and  gravitational  forces  in  the  wake  vorticity 
kinematics. 

PROPELLER  LOADING  VARIATION 

The  principal  parameter  varied  through  this  experiment  besides  die  propeller 
submergence  depth  was  the  propeller  loading  condition.  This  was  done  by  holding  the 
carriage  speed  constant  and  changing  die  rotational  speed  of  the  drive  shaft.  This  changed 
the  advance  ratio,  J,  representing  the  angle  of  attack  and  resultant  loading  on  the  propeller 
blades.  By  systematically  observing  the  wake  under  various  loading  conditions,  wake 
behavior  for  variations  in  propeller  jet  strength  could  be  evaluated. 

The  comparison  in  Fig.  8  (J  ranging  from  0.8  to  0.4)  shows  the  strong  dependence 
on  J  of  the  vortex  twisting  process.  The  most  obvious  effect  of  lowering  J  is  die  decrease 
in  streamwise  length  scale  (or  time  scale)  of  the  twisting  motion.  The  apparent  symmetry 
between  the  top  and  bottom  rows  of  vortices  indicates  negligible  influence  of  die  free 
surface  or  the  support  struts  at  this  deep  submergence  depth  (zWb/R=2).  The  influence  of 
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propeller  loading  on  the  vortex  interaction  process  was  therefore  pursued  through  closer 
examination  of  wake  hnages  under  various  J  conditions. 

Since  J  was  usually  varied  by  changing  propeller  shaft  angular  velocity  under  a 
constant  carriage  speed,  the  influence  of  unsteady  rotation  caused  by  the  universal  joints 
would  also  need  consideration.  The  influence  of  angular  velocity  on  J  can  be  determined 
from  the  definition;  J  =  U  /  (nD)  =  2jcU  /  (to,!)).  Computation  of  die  helical  vortex 
spacing  for  each  J  can  be  made  using  simple  geometry.  Estimates  of  die  propeller  blade 
load  can  also  be  made  using  the  open  water  propeller  loading  curves.  The  trends  in 
calculated  vortex  spacing  and  blade  loading  for  variations  in  J  are  shown  in  Fig.  9.  An 
estimate  of  the  tip  vortex  strength  can  be  made  using  the  blade  loading  by  assuming  that 
all  of  the  bound  vorticity  in  the  blade  is  shed  and  rolled  up  into  the  tip  vortex.  The  tip 
vortex  strength  is  known  then  from:  T  =  L/(pU)  (where:  L  =  lift  force  on  blade  /  span 
of  blade).  Note  that  the  lift  force  on  the  blade  can  be  estimated  from  the  open  water 
curves  for  a  known  J  value.  Combining  the  estimated  vortex  strength  and  separations  with 
the  Biot-Savart  law  gives  estimates  of  the  induced  velocity  of  the  vortices.  Computed 
trends  in  the  wake  vortex-induced  velocity  with  respect  to  J  are  plotted  in  Fig.  10. 

Induced  velocities  calculated  here  utilize  simplified  assumptions  of  low  vortex  curvature 
(essentially  2-dimensional)  but  illustrate  the  trend  of  a  full  3 -dimensional  analysis.  Note 
especially  the  increasing  trend  in  induced  velocity  between  die  vortices  as  J  is  reduced. 

The  higher  induced  velocities  at  low  J  values  speed  up  the  twisting  process  once 
perturbations  are  applied.  The  effect  is  similar  to  the  growth  in  numerical  instability  of 
filament  models  of  vortex  sheet  motions  as  the  filament  spacing  is  reduced  or  strength 
increased  as  described  in  Moore  (1971). 
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Using  this  induced  velocity  trend,  a  streamwise  length  scale  can  be  calculated  and 
compared  with  length  scale  measurements  from  the  video  images.  The  measurement  used 
for  comparison  is  the  streamwise  distance  from  the  propeller  plane  to  the  plane  where  a 
vortex  pair  has  first  reached  a  vertical  orientation.  Fig.  11  shows  a  comparison  of 
calculated  length  scale  (normalized  to  match  the  measured  data  at  J=0.5)  with  die  measured 
length  scale.  The  agreement  is  reasonably  good  considering  the  calculations  are  derived 
from  a  2-dimensional  model  of  the  vortex  interaction,  and  neglects  die  increase  in  induced 
velocity  caused  by  die  unsteady  shaft  rotation  at  lower  J  (see  Fig.  10)  and  decrease  in 
induced  velocity  caused  by  self-induced  velocities  of  curved  vortex  filaments.  The 
calculations  at  die  very  least  confirm  that  die  twisting  length  should  be  reduced  as  J  is 
reduced. 

An  additional  difference  in  die  vortex  interaction  was  noted  at  the  lower  J  values. 

At  J=0.4  the  larger  radial  gradient  in  streamwise  velocity  near  die  edge  of  the  propeller  jet 
causes  the  inboard  vortices  to  be  swept  downstream  past  their  original  partner  to  interact 
with  die  second  vortex  encountered.  Images  of  reduced  contrast  at  J=0.3  suggest  that  the 
vortices  are  swept  to  the  third  outboard  partner  downstream  before  interacting  to  swing  out 
of  the  jet  stream.  Turbulent  diffusion  of  the  dye  under  these  conditions,  however,  reduces 
the  detectability  of  organized  vorticity  in  the  high  velocity  fluid. 

It  should  be  noted  here  that  this  strong  wake  instability  should  be  greatly  reduced  in 
future  experiments  with  constant  propeller  rotation.  Despite  the  elimination  of  the  shaft 
velocity  perturbation,  future  experiments  are  expected  to  show  strong  dependence  of  the 
propeller  wake  evolution  on  J.  This  dependence  stems  from  the  inherent  instability  of  the 
helix  to  perturbations  (Widnall  1972),  and  the  decrease  in  vortex  spacing  as  J  is  decreased. 


In  particular,  it  will  be  interesting  to  note  the  variation  in  free  surface  influence  on  the 
breakdown  of  die  helical  form  of  the  propeller  wake  as  J  is  modified. 

PROPELLER  DEPTH  VARIATION 

Variation  in  die  submerged  depth  of  the  propeller  gave  rise  to  two  mutually 
interacting  effects  in  the  wake  of  die  propeller.  First,  variations  in  wake  vortical  structure 
downstream  are  influenced  by  the  presence  of  the  free  surface  boundary  condition.  In 
addition,  the  free  surface  elevation  is  modified  by  surface  piercing  model  support  struts  and 
local  flow  fluctuations.  These  local  flow  fluctuations  are  caused  by  proximity  of  either  the 
propeller  blades  or  their  tip  vortex  wakes.  The  two  depth  variation  influences  will  be 
treated  separately  in  die  sections  to  follow. 

Propeller  Wake  Influence 

The  influence  of  the  free  surface  on  the  vortical  structures  in  the  propeller  wake 
appeared  in  this  experiment  in  two  forms.  Variations  in  both  the  vortex  trajectory  and  die 
size/brightness  of  the  vortex  spots  were  noted  between  the  top  and  bottom  illuminations  of 
the  wake.  These  variations  were  amplified  as  J  or  die  non-dimensional  submergence  depth 
(z^JR)  were  decreased. 

A  study  of  Fig.  12,  showing  a  time  sequence  separated  by  0.1  seconds  at  J=0.5  and 
Zm/RsI.23  indicates  the  variation  in  vortex  trajectory  caused  by  proximity  to  the  flee 
surface.  In  this  sequence  the  vortex  labeled  "1”  moves  much  faster  over  the  top  of  its 
partner  labeled  "2"  than  the  corresponding  vortex  3  moves  below  vortex  4.  This  difference 
could  be  a  result  of  : 

1.  stronger  vorticity  in  vortex  1  or  2  than  in  vortex  3  or  4 
(possibly  due  to  reduced  dissipation) 
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2.  Stronger  shear  layer  on  top  with  respect  to  bottom  (possible 
slight  reverse  flow  on  top  due  to  zero  friction  boundary 
condition  at  surface) 

Note  that  the  movement  of  vortices  "2”  and  "4"  is  relatively  small. 

Also  notable  in  Fig.  12  is  a  distinct  difference  in  both  die  shape  and  brightness  of 
the  dye  spots  marking  the  tip  vortices.  As  mentioned  earlier  in  die  PROPELLER 
LOADING  VARIATION  section,  the  support  struts  show  no  influence  on  the  dye 
visualization.  The  authors  speculate,  therefore,  that  die  difference  in  dye  diffusion  between 
the  top  and  bottom  vortices  is  indicative  of  an  influence  of  the  free  surface  on  the 
turbulent  transition  of  the  vortices.  This  is  suggested  by  the  stronger  diffusion  of  dye  near 
die  edges  of  the  vortices  while  the  cotter  spot  stays  bright  longer.  Velocity  measurements 
in  future  experiments  will  shed  light  on  this  process. 

Free  Sal face  Influence 

As  mentioned  earlier,  die  free  surface  elevation  can  be  modified  by  the  presence  of 
several  types  of  disturbances. 

The  model  support  struts  piercing  die  free  surface  leave  waves  similar  to  those  left 
by  ships.  These  waves  are  dependent  on  the  Froude  number  based  on  the  speed  of  the 
carriage  and  chord  length  of  the  strut.  Observations  of  the  wave  patterns  behind  die  struts 
were  made  at  a  series  of  speeds  between  1  and  4  ft/sec  (0.305  and  1.22  m/sec)  at  the 
shallow  submergence  depth  of  z**/R=1.25.  During  these  observations  the  J  was  set  to 
1.045  corresponding  to  the  zero  thrust  condition  for  the  purpose  of  minimizing  propeller 
blade  and  wake  effects  on  the  free  surface.  These  effects  were  examined  separately  as 
described  in  the  next  paragraph.  The  waves  generated  by  the  struts  and  pod  are 
characterized  in  Table  1  below.  The  short  wavelength  transverse  waves  generated  at  1 
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A/sec  (0305  m/sec)  decreased  to  an  undetectable  amplitude  upstream  of  the  propeller  plane. 
In  addition  to  increased  dye  visibility  at  this  speed,  the  low  strut  influence  was  also  a 
factor  in  conducting  most  of  the  observations  at  this  speed. 

The  free  surface  effects  due  to  propeller  Made  and  wake  vortex  proximity  were 
amplified  expectedly  as  either  the  propeller  submergence  was  decreased  or  the  blade 
loading  was  increased.  At  2^JR=1.25  surface  depressions  caused  by  blade  passage  became 
deep  enough  (for  J<0.4)  to  draw  air  from  foe  surface  and  send  it  bubbling  through  die 
wake.  It  should  be  noted  that  this  ventilation  process  was  clearly  distinguishable  from  a 
cavitation  phenomenon  which  was  not  observed. 

Surface  seeding  proved  to  be  a  limitation  in  examining  die  deflections  at  the 
air/water  interface  caused  by  wake  vorticity.  Dust  particles  proved  to  be  very  helpful  in 
surface  visualization.  Strong  currents  created  near  the  surface  by  the  submerged  vorticity, 
however,  quickly  cleared  most  of  the  particles  deposited  there.  Remaining  particles  did 
indicate  correlated  free  surface  deflections  with  die  dye  marked  vorticity  below.  The 
dynamic  nature  of  the  video  image  is  necessary  however  to  fully  appreciate  this  correlation. 
Variations  in  the  laser  light  intensity,  caused  by  refraction  angle  changes  at  die  free  surface 
deflections,  were  also  sensitive  measurements  of  the  influence  of  submerged  vorticity  on 
the  free  surface  profile.  When  viewing  the  videotape,  varying  intensity  bands  in  the  laser 
light  sheet  can  be  traced  to  surface  deflections.  The  small  amplitude  of  these  surface 
deflections  under  the  conditions  of  this  experiment  are  usually  very  difficult  to  see 
individually  without  the  aid  of  these  light  rays. 
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WAKE  PERSISTENCE 


An  important  motivation  for  this  experiment  was  the  examination  of  die  persistence 
of  die  propeller  wake  a  long  time  after  the  propeller  had  passed.  Because  only  flow 
visualization  measurements  were  being  made,  and  the  luminescence  of  the  dye  decreased 
below  die  video  camera  threshold  after  approximately  40  seconds  at  low  J  and  60  seconds 
at  high  J,  no  long  term  wake  data  could  be  recorded  To  the  eye;  however,  die  fine  dye 
gradients  could  be  observed  for  over  4  minutes.  Comments  made  here,  will  therefore  be 
based  on  visual  observations  of  the  wake  made  by  the  authors  during  the  actual  testing. 
Two  characteristic  features  were  observed  in  the  far  wake  of  the  propeller  and  they  are 
described  below. 

The  first  feature  noticed  was  die  slow  migration  of  the  propeller  wake  jet  towards 
the  free  surface.  Once  wake  contact  was  made  with  the  free  surface,  higher  momentum 
fluid  quickly  moved  upward  forming  strong  axial  surface  currents.  These  currents  were 
visible  in  movement  of  both  the  dye-diluted  fluid  near  the  surface  and  dust  particles  on  the 
air/water  interface.  The  downstream  length  for  surface  contact  of  the  propeller  wake  with 
the  surface  (1,)  is  dependant  on  J  and  z**.  The  advance  ratio  essentially  determines  the 
wake  spreading  angle  p,  while  z**/R  defines  an  initial  depth  of  the  wake.  The  accurate 
measurement  of  1,  requires  a  precise  and  sensitive  criterion  for  evaluating  die  wake/free 
surface  interaction  point.  The  velocity  field  due  to  a  vortex,  however  prevents  the 
definition  of  some  range,  beyond  which  no  interaction  exists.  In  future  velocity 
measurement  experiments,  a  threshold  value  for  vortex  wake  induced  velocities  at  die  free 
surface  will  therefore  be  defined  to  signal  an  end  point  for  1,.  Measurement  of  P  was 
hampered  in  this  experiment  by  die  special  variation  in  the  wake  edge  caused  by  the 
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vortex  twisting.  A  more  analytic  evaluation  of  the  interdependence  of  these  variables  will 
therefore  require  additional  measurements  planned  in  die  future. 

Another  feature  noticed  in  the  far  wake  was  the  intermittent  occurrence  of  large 
turbulent  structures  long  after  the  propeller  had  passed  and  most  of  the  turbulent  motion 
had  diffused  to  low  intensity.  In  several  cases  these  structures  were  noticed  up  to  40 
propeller  diameters  downstream.  Velocity  measurements  may  show  additional  intermittent 
turbulent  structures  of  low  intensity  much  farther  downstream.  It  is  not  known  what  effect 
the  improvement  in  propeller  drive  smoothness  will  have  on  the  existence  or  strength  of 
these  persistent  turbulent  structures. 

DISCUSSION  AND  RECOMMENDATIONS 
The  experiment  described  here  included  video  recording  of  a  range  of  propeller 
loadings  and  carriage  speeds.  In  addition,  limited  video  coverage  of  the  3 -dimensional 
form  of  the  helix  was  also  obtained.  The  3 -dimensional  form  of  the  helix  (FK  6)  shows 
qualitatively  the  form  of  the  instability  present  in  the  data  reported  in  this  paper. 

Great  care  was  taken  to  align  the  universal  joints  in  setting  up  the  experiment.  No 
visual  sign  of  the  variation  could  be  seen  in  the  propeller  even  with  a  strobe  light. 

Detailed  measurements  described  in  the  HELIX  INSTABILITY  section  showed  however  a 
2.5  degree  amplitude  sinusoidal  variation  in  the  propeller  shaft  position  causing  blade 
loading  and  wake  helix  geometry  variations  of  significant  magnitude.  An  upgrade  in  the 
propeller  drive  design  utilizes  miter  gears  which  will  give  constant  velocity  at  the  propeller 
shaft.  A  retesting  of  the  propeller  with  the  new  drive  gears  will  be  carried  out. 

The  authors  suspect  that  once  the  shaft  velocity  variation  is  eliminated,  a  weaker 
instability  will  still  be  present  in  the  wake  helix  due  to  either  the  free  surface  effects  or 
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vortex  breakdown.  These  influences  should  be  more  prominent  in  future  experiments. 
Comparison  of  images  from  this  experiment  with  those  of  round  jet  shear  flows  suggest 
that  helical  vortex  geometry  in  the  presence  of  its  own  strong  shear  layer  may  be  unstable. 

The  implications  of  wake  helix  instability  could  have  a  significant  effect  on 
numerical  wake  analysis  codes.  The  generation  of  more  representative  input  mean  and  rms 
velocity  files  may  be  augmented  by  information  on  die  spatial  vorticity  distribution.  The 
presence  of  large  scale  vortical  structures  seen  far  downstream  in  the  wake  could  play  a 
significant  role  in  energy  storage  and  transfer  to  die  free  surface.  The  intermittent  nature 
of  these  structures,  however,  may  prevent  accurate  wake  modeling  by  mean  and  rms 
quantities  currently  being  used.  The  simultaneous  use  of  flow  visualization  and  Laser 
Doppler  Velocimetry  in  future  experiments  will  permit  cost  effective  evaluation  of  die 
importance  of  these  structures  and  their  interactive  nature  with  the  free  surface. 
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Fig  1.  Four  bladed  propeller  with  hollow  hub  and  dye  slotted  blades 


Fig  2.  Experimental  configuration 
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Video  frames  of  vortex  twisting,  J*0.5,  0.1  second  increments 
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Pig.  6  -  Three-dimensional  view  of  tip  vortex  1 
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Fig.  9.  Calculated  blade  loading  and  vortex  spacing  versus  J 
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Fig.  10.  Induced  velocity  variation 
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Table  1.  Free  surface  effects  due  to  struts  and  pod:  2^/11=1.25 


Carriage  Speed 
(ft/sec) 

1 


2 


3 


Wave  Pattern 

Short  1.5  in.  (4  cm)  wavelength  divergent  waves  due  to  the 
struts:  amplitude  <  0.5  in.  (1  cm). 

Diverging  and  transverse  waves  of  approximately  8  in. 

(20  cm)  wavelength  «  0.75  in.  (  2  cm)  amplitude. 

Dominant  transverse  wave  system  due  to  pod:  wavelength  * 
2  ft  (0.6  m)  and  amplitude  *»  0.8  in.  (2  cm)  due  to  pod. 
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Dominant  transverse  wave  system:  wavelength  «  3.5  ft 
(1.1  m)  and  amplitude  -  1  in.  (2.5  cm)  due  to  pod. 
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